A simple quantum model of spontaneous emission 
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We present a very simple model of a spontaneous emission from a two-level atom, interacting with 
a field of a finite number of states. Such a process is often said to occur because of the large number 
of equally- probable states of environment. We show that in our model increasing the number of field 
states may and may not cause a practically permanent emission, depending on the details of the 
model. We also describe how irreversibility emerges with growing number of states. Mathematical 
tools are reduced to a necessary minimum and hopefully can be well understood by undergraduate 
students. 



I. INTRODUCTION 



In the introductory course on classical thermodynam- 
ics, statistical thermodynamics is sometimes based on 
the concept of equally-probable microstatesi However, 
the microscopic mechanism which guarantees this equal 
probability is not discussed. 

Following this path, spontaneous emission may be said 
to occur because there exists large (or even infinite) num- 
ber of excited states of environment and only one excited 
state of the atom. It can be argued that all those states 
are equally probable and thus the equilibrium state is the 
atom in ground state and somehow excited environment. 
After presenting the model in Section UH in Section IIIII 
we present an example of a system, where the perma- 
nent emission is absent no matter how many states of 
the environment there are. 



On the other hand, in the typical academic course on 
quantum mechanics^, stationary problems constitute a 
major part. The only three places where dynamics ap- 
pear are the definition of time evolution operator, de- 
scription of two-state system and the time-dependent 
perturbation theory. The last is usually presented in a 
general way for arbitrary number of states. The exam- 
ple in which one could trace how the system properties 
change with varying number of states is missing. We give 
such an example in Section [IVI 



Thanks to the simplicity of the model, the approach is 
free of semi-classical arguments and does not involve ap- 
proximate methods like the perturbation theory. Thus, 
the model illustrates how fundamental principles of quan- 
tum dynamics works in a system more complex then a 
two-level one and can become a part of an introductory 
course on quantum mechanics, making it a little more 
oriented on dynamics than usual. Intuition build on this 
model may be of a great value when learning more ad- 
vanced topics, like open quantum systems or quantum 
field theory, where the huge or infinite number of states 
is an essential difficulty. 



II. BUILDING THE MODEL 

We assume that at the initial moment, the system con- 
tains single atom in excited state and the environment in 
the ground state. By spontaneous emission we under- 
stand the following situation: We prepare the atom in 
an excited state. We wait some time t w and measure if 
the atom is still excited. We repeat this procedure many 
times for similar, but not identical t w . We find that if t w 
is long enough, the atom is usually de- excited. Note, that 
this definition suits to the experimental reality. In par- 
ticular, t w is not the same in any two real measurements 
is obvious. 

We number all (orthogonal) basis states of the system 
by |0), |1), . . . , \N). We assume that N is finite. Thus, 
the space of states is an (N + l)-dimensional vector space 
and the Hamiltonian is an (N + 1) x (N + 1) Hermitian 
matrix (Hermiticity is required for the total probability 
to be conserved during time evolution). 

We choose |0) to denote the initial state (i.e. the state 
in which the atom is excited and the environment is in its 
ground state). The states \k) with k > correspond to 
different types of excitation of environment. These types 
of excitations we call " modes" , in reference to quantum 
optics. Notice, that all basis states are the states of the 
whole system, i.e. the atom and the environment. They 
differ in the part of the system, which is excited. 

Now we would like to write down the Hamiltonian of 
the system, which allows the system to flow from the 
state |0) to some state |fc). In quantum mechanics time 
evolution of the system is described by the time-evolution 
operator Ut, 



(i) 



With the spectral theorem, we immediately see that in 
the eigenbasis of H the operator Ut is a diagonal matrix 
with matrix elements 



(e n \U t \e r , 



a iE n t/h 



(2) 



where n,n' £ {0,1, ... ,N}, |e„) are eigenstates of H, 
E n are corresponding eigenvalues and S denotes the Kro- 
necker delta. In arbitrary basis for small At we have 



U At w I + iHAt/h 



(3) 
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and we see that if the Hamiltonian matrix element 
(fc|-ff|0) is non-zero, \ip(t + At)) has different component 
\k) then \ip(t)). Otherwise, (k\if>(t + At)) = (k\ip(t)). 
This is valid for all t, so (fc|iJ|0) must be non-zero if 
evolution from |0) directly to \k) is possible. We denote 
(k\H\0) — ak- Due to Hermiticity, also (0\H\k) = a* k 
must be non-zero, so the reverse process is also possible. 
This is true for all k > 0. 

We assume that different modes are independent, i.e. 
the Hamiltonian does not mix them, (k\H\k') = for 
k ^ k! ' . This assumption is not necessary, but simplifies 
the model. We denote (k\H\k) (k £ {0, 1, ... , N}) by e k . 

The probability that the atom is excited after time t is 
P(t) = |(0|V>(*))| 2 - Using Eq. (JT) and ©, we can write 
it as 



P(t) 



N 

E 

n=0 



K0|e n )| 



2 iE n t/H 



(4) 



III. DESTROYING STEREOTYPES 

Now we would like to describe with our model two 
systems, where the spontaneous emission is not present. 
They are counterexamples for some widespread stereo- 
types. 

The first example is the case N = 1 (there is only one 
excited state of environment), E\ = Eq (mean energies of 
the state with excited atom and of the state with excited 
environment are the same) and arbitrary ct\ — a. Then 
H is 2 x 2 matrix, 

Hi = £0 (|0)<0| + |l>(l|) + a |0)<l| + a*|l>(0| (5) 

e a 
a* e 

It is easy to compute eigenvalues and eigenstates from 
characteristic equation. Then from Eq. (j4]) follows 



(6) 



Pi(t) = cos 2 ( l -ft 



This mean that no spontaneous emission is present. 
The atom oscillates between the ground and the excited 
states. We can say it is at resonance with the field. Stu- 
dents can individually check that similar result holds for 
arbitrary E\ (i.e. out of resonance), but the amplitude 
of oscillation is smaller and the atom is never certainly 
de-excited. 

This example is a special case of the familiar 2-level 
system. It destroys the stereotype, that "going to lower 
energy" is prescribed in any quantum theory. 

As the second example we consider the case of N iden- 
tical states, i.e. Ek — Eq and ak = a, for some given N. 
The Hamiltonian now has the form 



N 



H 2 = e £l fc X fc l+EH X fc l+ a *l fc X°l)- (7) 



This example can be solved analytically by a simple 
change of the basis. We denote 



1 N 

i) = ^Ei fc ) 



(8) 



The new basis consists of |0) = |0), |1), and any orthog- 
onal to them (and each other) and normalized states 
1 2), . . . , \N). The part of Hamiltonian proportional to 
the identity matrix X^fcLo I^X^I l°°ks the same in any 
basis. Using Eq. © and the facts that (0\k) — is the 
condition of orthogonality to |0) and J2 k=1 (k\k) = is 
the condition of orthogonality to |1), the following part 
can be expressed as 



N 



N N 



^2a\0)(k\+h.c. = ^2^2a\0){k\k)(k\+h.c, (9) 

fe=l k=l k=0 

= y/Na\0)(l\+h.c. 

(h.c. denotes "Hermitian conjugate"). Thus all states or- 
thogonal to |0) and |1) are eigenstates with eigenvalue £q. 
This means that we only have to find two last eigenstates, 
i.e. diagonalize remaining 2x2 submatrix. Further cal- 
culation and the resulting evolution is similar to the case 
N = 1. The probability that the atom is excited is 



P 2 (t) = cos 2 



N\a\ 



(10) 



No spontaneous emission is present, no matter how large 
N is. This example destroys the stereotype that mi- 
crostates are always equally probable and the huge num- 
ber of possible states affirm that spontaneous emission is 
present in the system. 



IV. CONDITIONS FOR SPONTANEOUS 
EMISSION 

To understand when spontaneous emission occurs we 
analyse Eq. ((3]). We notice that the expression for P(t) 
resembles a bit a square of module of some Fourier series. 
Consider the special case when TV = 2M. To parametrize 
eigenvalues and eigenstates we use m £ {— M, . . . , M} in- 
stead of n £ {0, . . . 2M}. We assume that the eigenvalues 
are equally-spaced in the section eq ± D, 



E,, 



so 



m 



D. 



and that 



(11) 



(12) 



fc=0 



fc=i 



We would like to stress, that now we make assumptions 
about the eigenvalues of the whole system and we have 
not yet discussed the Hamiltonian matrix elements Ek 
and ak in our initial basis. 
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From Eq. §4§ we have 



M 



,ie t/H 



E l<oi 



-M 
M 



(13) 



l(0|e o )| 2 +E 2 l 



e m cos 



Eq. (fT3"|) states that d P-$(t) can be expand in a (finite) 
Fourier series. Thus it is a periodic function with period 



M M 
T = 2ttH— = h—. 
D D 



(14) 
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However, this does not mean that spontaneous emission 
is not present in the system. First, in reality the dimen- 
sion of space of states is surely huge, possibly infinite. 
For large M, T is also large. There exists such a big M 
that T is too long to be ever measu red. What is more, 
even when T is quite small, \J Pit) may be measurably 
non-zero only for a very short time r in every period. If 
unavoidable differences in t w in different measurements 
are significantly larger then r, spontaneous emission is 
observed. 

Note that the coefficients of Fourier expansion are de- 
pendent only on scalar products of initial state |0) with 
eigenstates |e m ). The only properties that they must ful- 
fil for all choices of £k and ak are that they are positive 
real numbers summing to the unity. To prove this we 
show the algorithm which finds and 

Choose |(0|e m )| 2 > 0. E m fulfils Eq. (HU). In its 
eigenbasis, H is a diagonal matrix with elements 

Em- 

|0) may be chosen to have real, positive compo- 
nents in the eigenbasis of H . Their moduli arc de- 
termined by the choice of |(0|e m )| 2 . 

Orthonormalize {|e_A/), . . . , |e_i), |0), |ei), . . . , |ejw)} 
using Gram-Schmidt procedure, starting from |0). 
Find H in this basis. 

Consider the subspace orthogonal to |0). Denote H 
reduced to this subspace by H re d- Find the eigen- 
basis of H re d- 

5. Change the basis such, that |0) stays invariant and 
the following basis states are eigenstates of H re d 
extended to the original space. Calculate H in this 
basis. 

6. The result is the Hamiltonian matrix which has the 
desired form. 

Realizing this algorithm numerically for some concrete 
|(0|e m )| {e.g. in the example considered below) may be 
an instructive exercise for students. 

The corollary is that there exist very wide range of 
possible P(t). For example, we can choose 



4. 



l(0|e m )| 5 



1 



for any m. This is not a blind guess. We know, that 
the Fourier transform of the Dirac delta function is con- 
stant, so we expect that with Eq. (fT5|) y/ P3 (t) will be 
considerably different from zero only for t being close to 
integer multiple of T. The results are plotted in Figure 
1. It can be seen, how period becomes larger when M 
increases. Simultaneously, the time after which the emis- 
sion happens is approximately independent of M. Thus 
for big enough (but still finite) M, the probability that 
after long and not precisely determined time t w the atom 
is excited is not measurable. 




2M + 1 



(15) 



Figure 1: Probability P{i) that the atom is in excited state 
after time t when the relation of the initial state to eigenstates 
is given by Eq. (|15p and the energy eigenvalues are given by 
Eq. (fTTj) . Total number of states equals 2M + 1. 



V. CONCLUSIONS OF THE MODEL 

We propose a simple quantum model for the system 
consisted of the two-level atom and ./V-modal field. We 
have shown, that large number of modes is a necessary 
but (contrary to widespread opinion) not sufficient con- 
dition for occurrence of spontaneous emission. Whether 
the emission is present in the system or not, depends on 
the interaction with the field and the modes mean ener- 
gies, which determine the coefficients |(0|e m )| 2 in Eq. j4|. 

We have shown the algorithm of building a Hamil- 
tonian which cause P(t) to be desired function if only 
square root of this function has a finite expansion in the 
cosine Fourier series (with any period). This is an ex- 
tremely broad set of functions. This shows that even a 
simple two-level atom can undergo various types of an 
evolution when is placed in a complex environment. 
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